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Abstract 

In this paper we investigate backreaction of excitations on a planar domain 
wall in a real scalar field model. The backreaction is investigated in the cases 
of homogeneous, plane wave and wave packet type excitations. We find that 
the excited domain wall radiates. The method of calculating backreaction 
for the general forms of excitations is also presented. 



1 Introduction 



The presence of topological defects in the field theoretical models with a 
degenerate vacuum is an important aspect of the structure of these models. 
Dynamics of topological defects is to be extracted from nonlinear equations 
describing evolution of fields they are composed of. The results, eventhough 
it is a formidable task to get them, are of great interest for particle physics 
(e.g. dynamics of a flux-tube in QCD [1]), for field theoretical cosmology (e.g. 
cosmic strings [2], [3]) and for condensed matter physics (e.g. domain walls 
in magnetics, vortices in superconductors or in superliquid helium, defects in 
hquid crystals [4], [5], [6], [7], [8]). 

This paper is devoted to dynamics of domain walls governed by a Poincare 
invariant wave equation. Domain walls appear in the models with the non- 
trivial zeroth homotopy group of the vacuum manifold. In the papers [9], 
[10], [11] two main approaches to the dynamics of domain walls have been 
presented. The first one is the polynomial approximation, the second one is 
the expansion in the width of the wall. Calculations made in the framework 
of these two methods indicate existence of an oscillating component in the 
width of the domain wall. This suggests that the domain wall can radiate. 
The radiation was also observed in computer simulations [12]. 

The problem we adress ourselves to is connected with this special aspect of 
dynamics of domain walls, namely the radiation. Wc calculate the radiation 
with the help taken from [13], [14] where analogous was considered in the 
abelian Higgs model. The method consists of two steps. The first one is to 
find the excitations of a static domain wall. They are investigated in the 
linear approximation and treated as small corrections to the basic domain 
wall field which is localized on the wall. The second step consists of looking 
for the effects of this excitations on the evolution of the domain wall, i.e. a 
backreaction. This procedure can be treated as the expansion of the domain 
wall field in the amplitude of the excitation A (A<Sl). The zeroth order term 
is then the static, exact planar domain wall solution (j)Q, the first order term 
is the excitation 0i and the second order term is the backreaction 

(l)^(j)o + A(j)i + A^(f)2 + 0{A^). (1) 

The terms of the higher order can be interpreted as more complicated effects, 
e.g. the third order term as the selfinteraction of the excitation. We shall 
discuss three main cases of the excitations, namely the most simple of them 
- the homogenous excitation, the second case - a plane wave type excitation 
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and the last one - a localized excitation. In this paper we shall consider the 
dynamics of the domain wall in the simplest scalar field model which has a 
potential with two degenerate minima. 

The plan of our paper is the following. In the next Section we present the 
calculation of the backreaction of the homogeneously excited domain wall. In 
Section 3 we describe the method of calculating backreaction for the general 
forms of excitations. Section 4 is devoted to the detailed analysis of the 
backreaction in the cases of plane wave and wave packet type excitations. In 
Section 5 we summarize the main points of our work. 

2 Radiation from a homogeneously excited 
domain wall 

We consider the planar domain wall in the model defined by the Lagrangian 
density, 

/: = ~rj,ud''^d''^-^{^'-vY: (2) 

where r/^i, =diag(— 1, 1,1,1) and A, u are positive constants. The correspond- 
ing evolution equation for the scalar field $ has the form: 

d^'df,^ - 2A($2 _ v^)^ = 0. (3) 

Exact, static solution representing the domain wall localised around the 
{x^, a;^)-plane is given by the formula: 

4>o = V tsmh^ax^) , (4) 

where a — y/Xv^. 

Let us rescale the scalar field $ and the space-time coordinates x^, 

x^ — ax^, 

where (p{x^) and are dimensionless. Moreover, for x° and x^ we shall use 
the following notation. 



(5) 



In the rescaled variables the evolution equation and the static solution take 
the form: 

d'^di.cl) - 2(02 - 1)0 = 0, 

(7) 

00 = tanh(0, 

where = The first step of our considerations is to find excitations of 
the planar domain wall. As it was stated in the Introduction the excitations 
are considered as small corrections to the basic static domain wall solution. 
Thus the corresponding equation for the excitation field 0i we shall obtain 
as the linear approximation to the initial evolution equation. Inserting the 
expansion (1) into Eq.(7) and keeping the terms of the first order in the 
expansion parameter A we obtain: 

d%(l), - 2(30o' - 1)01 = 0. (8) 

The planar domain wall distinguishes the direction perpendicular to the wall 
plane, given in our case by the coordinate lines of ^. Small correction 0i 
do not change this asymmetry. Therefore we assume a special form of the 
perturbation in which dependence on the coordinate ^ is separeted, 

<t>i^mx{r,x\x'). (9) 

Inserting the above Ansatz into Eq. (8) we obtain the system of two equations 
for the functions and x with a separation constant c, 

d^daX = cx, 

where a=0,l,2. Prom that follows a restiction for constant c. If c were less 
than zero the solution of the second equation would increase with the time 
infinitely and so wouldn't meet the main requirement of the perturbative 
calculation. 

Let us solve the first equation of the system (10). It can be transformed 
into the following form: 

1 d'^i^ 3 , ,c ^, , , , 
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It is the special case of Schrodinger equation with generahzcd Poschel- Teller 
potential. Solutions of this equation can be found in [15]. There exist two 
bound states enumerated by n=0,l, 

n = V^o(0 = 3SiN, Co = 0, 

(12) 

^ = 1 MO-^,, c, = 3. 

Both are localized on the wall in the sense that they exponentially decrease 
for large | ^ |. The first, even solution ■00 (0 can be interpreted as a small 
displacement of the domain wall, 

= ^o + V'oX = tanh^H - tanh(4 + x)- (13) 

cosh 4 

This is the zero mode related to the translational symmetry of the model. 
The excitation given by the second, odd function il^i does not posses such 
interpretation. This solution we shall accept as the proper bound state of 
the domain wall. 

The second equation of the system (10) is of the wave type. In this Section 
we shall consider the special case of excitation, homogeneous on the whole 
wall and given by the solution independent of the coordinates x^,x'^, 

X{r) = Acos{y/3T + 5), (14) 

where 5 is a constant phase which we shall put equal to zero, and A is a 
constant amplitude. As mentioned above the described procedure assumes 
that A is sufficiently small. Finally, the field of the homogeneously excited 
planar domain wall is given by the formula: 

= tanh^ + A 2^cos(y3r). (15) 

cosh 4 

Prom this formula one can see that the excitation introduces periodic changes 
of the thickness of the domain wall. 

The second step in our procedure is to find the backreaction. Inserting 
the expansion (1) into the Eq.(7) and keeping all terms of order we obtain 
the equation for the backreaction 02, 

d^d^<j>2 - 2(30o' - 1)02 = 6020i'. (16) 
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It is an inhomogeneous equation with the r.h.s. including the square of the 
excitation component. We will denote the inhomogeneous term by = 
6^001 • For a homogeneous excitation we are dealing with, is a function 
of the two variables (^,t). Inserting the formulae for the functions 0o and 
(pi we obtain: 

„ . sinh^f ,9sinh^f , /- . , . 

N^Ni + N2^SA2 — r^ + SA^ — -^cos(2V3t). (17) 
cosh 4 cosh ^ 

We shall solve Eq.(16) in two steps, considering each of the two parts of the 
inhomogeneous term separately. At the first step we consider the equation: 

'^<P'(i>2 / ,9.. X , 3 .oSinh^^ , , 

-^-(3ta.i=5-lM. = -^^;^. (18) 

The general solution of this equation can be found by the standard Green's 
function technique, see e.g. [16]. As the two linearly independent solutions 
of the homogeneous part of Eq.(16) we take: 

(19) 

/,(0 = |sinh(20 + |tanh^+|^. 
As the Green's function we take 

G{^,x) = 2h{x)h{mi - ^) - Vi{0f2{x)m - ^) - 0{-x)]. (20) 

The Green's function was choosen in such a manner as to obey the condition 

G'(^ = 0,x) = f\xeR. (21) 

Such a choice ensures that an inhomogenuity won't produce any displacement 
of the domain wall as a whole. The general solution of Eq.(18) has the form 

/+0O 
G{C,x)h{x)dx, (22) 
-oo 

where 



h{x) 



3 sinh^ X 



2 cosh X 



5 ^' 
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and a, b are constants. Formula (22) gives: 



MO = A'iciOMO + C2(0/2(0], (23) 
ci(0 = - 2 /o^ f2{x)h{x)dx, , . 

The function fi is even while /2 is odd. Because f^iO exponentially grows for 
^ ±oo the coefficient fTinction C2 has to vanish in this limit. Therefore we 
have to put 6 = while a is still arbitrary. But we put a to be zero too because 
keeping it nonzero would amount to including uninteresting solution of the 
homogeneous equation. In the second step wc have to solve the Eq.(16) with 
the second inhomogenous term N2 containing a periodic time-dependence. 

The perturbatively obtained Eq.(16) for the beackreaction has radiation 
type solutions which do not vanisch for large ^ as will be shown further. In 
this case we adopt the Helmholtz condition which states that for | ^ |^ 00 
only outgoing radiation waves are present. 

Hence we consider the following form of the solution: 

(j>2 = exp (-z2V3r) + exp (z2x/3r)], 

where the positive and negative frequency components are related by complex 
conjugation, 

This Ansatz leads to the following equation for the functions (f±, 



1 d 



■2 



+ (4 + 



2de cosh^^ 



3 sinh^ ^ 



which we solve analogously as Eq.(18). As two linearly independent solutions 
we take 

g,{0 = cos(2y20 - V2tanhesin(2v^O + I'^f, 

(26) 

92(0 = sin(2V20 + x/2tanhecos(2v^e) + I'^^f- 
The Green's function is given by the formula 

GiC, x) = -l=[^2(0^i(^)^(^ - ^) - 9i{092{xm^ -x)- 9{-x))]. (27) 
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The general solution has the form 

^±{0 = [dliO + «±]^i(0 + [d2iO + b^MO, (28) 

where ^ 

MO ^--^l 92{x)h{x)dx, (29) 

d2{0 = ^ £^ gi{x)h{x)dx, (30) 

and a±, b± are constants. 

In order to satisfy the Helmholtz condition imposed on the solution, we 
have to find solutions with the appriopriate asymptotics given by the formula 
below 

(^±({-> ±oo) ~exp[±i2y2 I ^ I]. (31) 

We consider separately the regions ^ > and ^ < and choose the solutions 
which will have the proper asymptotic behaviour given by the formula (31). 
Next, we impose the matching conditions at the point ^ = 0, i.e. the conti- 
nuity conditions of the solution and its first derivative. The solutions found 
in this way are given by the formula 

V± = d^{Ogl{0 + [d2{0 ± tdi{oo)MO, (32) 

where 

di{oo) = lim di{^) 
and is finite. Finally, the asymptotic form of the solution is the following 

(p2{^ ±oo) ~ ±V3di{oo)A'^ cos(2^ | ^ | -2^r ± /?), (33) 

where P = arctan \/2, G (0, |). 

Let us summarize results of our calculations. The full solution consisting 
of the static domain wall, the excitation and backreaction is the following, 

4> = tanh^ + cos(V3r) + A^[c,{Oh{0 + ^2(0/2(01+ 

(34) 

+ lA^[ip+{^) exp (-i2x/3r) + ^_{^)exp{i2y/3T)]. 

Backreaction consists of two terms. The first one, independent of time, con- 
sists of the static reaction of the domain wall to the homogeneous excitation 
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while the second one depends on time. Its asymptotics is given by the for- 
mula (33) and describes the radiation with the energy fluxes given by the 
Poynting vectors: 

S± = 6V3dl{oo)A^v'^a'^ sin^k^x^ ± (3)k±, (35) 

where 

A;^ = (2^3, 0,0, ±2^2). (36) 
and the signs ± correspond to the limits +oo and —oo respectively. 

3 The method for the general forms of exci- 
tations 

Our goal in this Section is to present the method of calculating backreaction 
for a general form of excitation (pi defined by the formula: 

(t>i^mxir,x\x'), (37) 

where 

sinh^ 

HO = — 727' 
cosh 4 

and the function x is any bounded solution of the wave equation 

d-daX = 3x. (38) 
Let us recall the equation for the backreaction, 

d%(t>2 - 2(3tanh' ^ - 1)02 = 6(^o0i'- (39) 
The first step is to define the operator by the formula: 

Then the Eq.(39) takes the form: 

d'^dah + L^h - 402 = 60o0i'. (41) 
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The second step is to solve the eigenvalue problem for this operator. We 
shall find all solutions of the equation: 

L^MO = A^a(0, (42) 

where A is the eigenvalue and '0a(O - the normalized eigenfunction corre- 
sponding to this eigenvalue. The backreaction 02 we treat as the expansion 
in the eigenfunctions of given by the formula: 

= JaA(£")V'A(0, (43) 

where the coefficients a\ depend on the coordinates x". Inserting this expan- 
sion into the Eq.(41) we obtain 

^, ^A'(0[^"4 + A' - AMx^) = NiC, x"), (44) 

where iV(^,x°) = 60o0i^- Multiplying the equation above by the eigenfunc- 
tion if^xi^) and integrating over the full range of variation of ^ and using the 
ortonormality condition for the system of the eigenfunctions '0a (0 we obtain 
the system of differential equations for the coeffcient functions ax{x"-), 

(a«4 + A-4)aA(F) = /iA(F), (45) 

where 

/+0O 
-oo 

The last step in this procedure is to solve Eq.(45). It may be done by 
the standard method of the Green's function. Let us denote by Gx{x"') the 
retarded Green's function of the operator, 

^A = d^da + A - 4. (46) 

Then the solution of the inhomogeneous equations (45) is given by the for- 
mula 

axix") = / Gxix" - i''")hx{i''')di'\ (47) 

In this solution we have dropped the homogeneous part because it does not 
contain any information about the excitation. 
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The results of the calculations are the following. The spectrum of the 
operator consists of two parts. The discrete part, for which A > 0, contains 
two eigenf unctions, 

The continuous part, for which A < 0, inchidcs twice degenerate subspaces 
corresponding to the eigenvalues A = —k^{k G R+) spanned by the eigen- 
functions: 



V^,(^) (0 = nik) [{k' - 2) cos(A;0 - 3k sm{kO tanh ^ + '-^] 



cosh^ 



(^) ^ n{k) [{e - 2) sm{kO + 3k cos{ki) tanh i + 

where 



(49) 



n{k) = [7rA;(A;2 + l)(fc2 + 4)]' 



Therefore we have to find the retarded Green's function for the following 
operators: 

Di=~d'^~da (A = 4). 

D2 = 9"4-3 (A = l), (50) 

= - - 4 (A = -A;2). 

They are obtained by the Fourier transform method and given by the formu- 
lae: 

^1(5") = 



G,{~z^) = r dKK^^^^^^MKz), (51) 

The procedure presented above enables us to find the backreaction for the 
general form of excitation of domain wall. We shall use it in the next Section. 



where 
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4 The backreaction for the plane wave and 
wave packet excitations 

In this section we analyse the backreaction of the excitation of the plane 
wave and wave packet type. In the first case x is given by formula: 

X = Acos{uj{kl,kl)T - klx^ - klx^), (52) 

where uj{kl kl) = ^{klf + {klY + 3. 

In the second case we consider the approximate solution of Eq.(38): 

X = Aexpi- ^"^ \ sin(V3r) = Aw{x) sin(v^r), (53) 

where A and A are constants. We assume that A ^ 1 and that the wave 
packet has momentum cutoff k ~ A^^. Then we may neglect for the finite 
time interval < r <~ A the spreading out of the wave packet, which is of 
course present in exact wave packet solutions of Eq(38). 

l.The plane wave case. 
The inhomogeneous terms in the backreaction equation has the following 
form: 

iV(^) (r, ~x\ x\ i) = 3^2— -^[1 - cos(2a;(A;^, kl)T - 2k',x' - 2klx')] (54) 

cosh t, 

It is convenient to pass to the Fourier transform with respect to the coordi- 

#(e,x^^^r) = 



nates x^,x^, r: 



Jdk^ Jdk^ Jdcue^p [-iiur + ik^x^ + ik'^x^U^^i^, k\ k^,u), 

(27r)5 

N^^\i,x\x^,T) = 
^ Jdk^ Jdk^ Jdcue^p [-iojT + ik^x^ + ik'^x^\N^^\i, k\e,u}). 



(55) 



(56) 



In this case, 

iV«(e, k\ k\ a;) = 3(27r)M2^[5(u;)5(A;i)^(fc2) 

-\5{uj + 2Lo{kl kl))5{k^ + 2kl)5{e + 2kl) (57) 
-\5{u - 2u{kl, kl))5{k^ - 2kl)S{k^ - 2kl)]. 
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The equation for the Fourier transform 02^^ has the following form: 

- + ^ - 2(3 tanh^ ^ - l)]^^'^^, k\ k\ a;) = N^'\^, k\ k\ uj). (58) 
It is clear from the above equation that the solution has the form: 



{i,k\k\u)^A\27^n^^\i)5{u)5{k^)5{k^) 



^"-l^OSiuj + 2uj{kl kl))5{k^ + 2kl)5{k'' + 2kl) 



(59) 



ip'l\m^ - Mkl kimk' - 2kl)5{k^ - 2kl)]. 



The negative and positive frequency components are related by complex con- 
jugation, 

if^l^ = (60) 



while (/Pg^^ is real valued. The functions fg^^, ip± have to satisfy the following 
equations: 



1 (f 3 

+ 



2de ' cosh^^ 
1 (P 3 



- 2 



+ 4 



(1) 3 sinh^ ^ 



5 



2 cosh^ ^ 



(1) 3 sinh ^ 
" 2cosh^e' 



(61) 
(62) 



2 cosh^ i 

The function (/^g^"* corresponding to the frequency = contains the informa- 
tion about the static backreaction of the domain wall while two remaining 
functions describe the dynamic backreaction. We shall concentrate on these 
two functions only. They must satisfy identical conditions as the functions 
y9± in Section 2 and they also satisfy the same equation, see formulae (25 - 
31). Therefore, the solutions arc given by the formula (32). The asymptotic 
form of the radiation part of the backreaction is then the following: 



#(x^x^e,r)~ 

^^/2,A^dl{oo) cos [2/c^x^ + 2klx^ + 2^/2 \ ^ \ -2uj{kl, k^r ± /?], 



(63) 



where /3 = arctan -\/2, /3 e (0, |). The corresponding wave vectors have the 
components: 

(k^) = {2u;{klkl),2kl2kl±2V2), (64) 
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where the signs ± correspond to the hmit in the ±00 respectively. The energy 
fluxes due to these waves are given by the Poynting vectors: 

S± = 6dl{oo)A%^a'^uj{kl, kl) sin^ik'ix^ ± f3)k±. (65) 

2. The wave packet case. 
Calculations of backreaction are carried out in the analogous steps as in plane 
wave case. Instead of formulae (54), (57) we have now: 

7V(2) (e, x\ x\ t) ^ ZA'^^^^y{x)[l - cos(2V3r)], (66) 

cosh ^ 

A^(2)(^,A;^fc^a;) =3A^— -4«;(fc)[5(a;)--<5(a;+2V3)--%-2V^)], (67) 

cosh 422 

where ^ 

w{k) = — Jdx^ Jdx^ eKp[—ik^x^ — ik'^x^]'w^{x), 

A; = ^(A;l)2 + (/c2)2_ 

We next obtain the same equation for the Fourier transform 02^^ (^, k^, k'^,u}) 
of #(e,xSx2,T): 



^'-^' + :£-2(3tanh='^-l) 



. 0?)=iV(^). (68) 
Analogously to the previous case the solution can be written as 
= A''^^^wik)[^^^\0Siu;)-l^^^\0S{u;+2V3)-l^^^^ 

(69) 

(2) (2) 

In the exact solution the functions if g ,Lp-^ dependent on k, but in the case 
at hand the expression in square brackets on the left hand side of Eq.(68) 
can be simplified. Namely, we may neglect the term k"^ because of the cutoff 
A » 1. This simplification leads to the same set of equations (61), (62) for 

(2) (2) 

the functions ipQ ,ipj. as in the previous case and the solution given by the 
formula (69) is then the approximate one. The asymptotic conditions (31) 
remain unchanged. Thus the asymptotic form of the radiation part of the 
backreaction is the following: 

0^^^(C ^ ±oo,x^^^T) ~ ^V3di(oo)A'^cos(2V2 \ ^ \ -2VSt ± (3)w'^(x). 

(70) 
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The corresponding wave vector has the components: 



A;^ = (2\/3,0,0,±2\/2). 



(71) 



The energy flux is given by the formula: 



(72) 



5 Remarks 

1. Let us summarize the main points of our work. We have presented the 
calculations of the backreaction in the cases of the homogeneous, plane wave 
and wave packet type excitations of the domain wall. We also have described 
the method enabling us to analyse the more general cases of the excitations. 
The main result of our work is the existence of the long range component in 
the backrection which is interpreted as the radiation from the excited domain 
wall. The frequency of the radiation, given by 02, (cf. formulae (63), (70)), 
is twice of that of the excitation function (cf- formulae (52), (53)). 

2. The idea of the expansion in the amplitude of excitation, which all the 
calculations were based on, was apphed to the simplest model of the real 
scalar field and the potential V{(f)) — |(0^ — f^)^. It seems that without 
much trouble this method could be applied also to the other field-theoretical 
models containing domain wall configuration. On the other hand, we should 
remember that this method is based on the linear approximation, what im- 
phes that it is rehable for small amphtudes of the excitations only. In order 
to consider stronger excitations we have to work out another method which 
will take into account the nonlinearity of the evolution equation in a better 
manner. One could for instance use the polynomial approximation in the 
vicinity of the domain wall and the proper asymptotics at the infinity and 
smoothly match them in the intermediate region. 
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